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We present the first simultaneous measurement of the ratio of branching fractions, R — B(t — > 
Wb)/B(t — ► Wq), with q being a d, s, or b quark, and the top quark pair production cross section o t i 
in the lepton plus jets channel using 0.9 fb _1 of pp collision data at ^/s = 1.96 TeV collected with 
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the DO detector. We extract R and a t j by analyzing samples of events with 0, 1 and > 2 identified 
b jets. We measure R = 0.97l°;os (stat+syst) and a t i = 8.18±°;8° (stat+syst) ± 0.50 (lumi) pb, in 
agreement with the standard model prediction. 

PACS numbers: 13.85.Lg, 13.85.Qk, 14.65.Ha 



Within the standard model (SM) the top quark decays 
to a W boson and a down- type quark q (q — d, s,b) with 
a rate proportional to the squared Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element, \V tq \ 2 Under the 
assumption of three fermion families and a unitary 3x3 
CKM matrix, the | Vt q | elements are severely constrained: 
\Vtd\ = (7.4 ± 0.8) • 10" 3 , \V ts \ = (40.6 ± 2.7) • 10~ 3 and 
\V tb \ = 0.999100tgXiM [2]. However, in several exten- 
sions of the SM the 3x3 CKM submatrix would not 
appear unitary and \V tq \ elements can significantly devi- 
ate from their SM values. This would affect the rate for 
single top quarks production via the electroweak inter- 
action [3| and the ratio R of the top quark branching 
fractions, which can be expressed in terms of the CKM 
matrix elements as 

R _ B(t^ Wb) = I Vtb I 2 

B(t ^Wq) | Vtt | 2 + | Vts | 2 + | V td | 2 ■ 

A precise measurement of R is therefore a necessary in- 
gredient for performing direct measurements, free of as- 
sumptions about the number of quark families or the uni- 
tarity of the CKM matrix, of the \Vt q \ elements via the 
combination with future measurements of the single top 
quark production in s and t channels [4| . 

In this Letter, we report the first simultaneous mea- 
surement of R and the top quark pair (tt) production 
cross section a^. R was measured by the CDF and DO 
collaborations [5, 6]. The simultaneous measurement of 
R and a t i, in contrast to previous measurements , al- 
lows one to extract a t j without assuming B(t — » Wb) = 1, 
and to achieve a higher precision on both quantities by 
exploiting their different sensitivity to systematic uncer- 
tainties. 

The current measurement is based on data collected 
with the DO detector Q between August 2002 and De- 
cember 2005 at the Fermilab Tevatron pp collider at 
y/s = 1.96 TeV, corresponding to an integrated luminos- 
ity of about 0.9 fb _1 , approximately four times larger 
than that of our previous measurement The analysis 
uses the top quark pair decay channel tt — > W + qW~q, 
with the subsequent decay of one W boson into two 
quarks, and the other one into an electron or muon and 
a neutrino. This is referred to as the lepton plus jets 
(^+jets) channel. We select a data sample enriched in tt 
events by requiring > 3 jets with transverse momentum 
Pt > 20 GeV and pseudorapidity \r)\ < 2.5 one iso- 
lated electron with pt > 20 GeV and \r]\ < 1.1 or muon 
with pr > 20 GeV and |r7| < 2.0, and missing transverse 
energy $ T > 20 GeV in the e+jets and $ T > 25 GeV in 
the /x+jets channel. Additionally, the leading jet pr is 



required to exceed 40 GeV. Events containing a second 
isolated electron or muon with pr > 15 GeV are rejected. 
The lepton isolation criteria are based on calorimeter and 
tracking information. Details of lepton, jets and^T iden- 
tification are described elsewhere [lOj . 

We identify 6-jets using a neural-network tagging al- 
gorithm [ll[ . It combines variables that characterize the 
presence and properties of secondary vertices and tracks 
with high impact parameter inside the jet. In the simu- 
lation, we assign a probability for each jet to be 6-tagged 
based on its flavor, px, and r/. These probabilities are 
determined from data control samples, and can be com- 
bined to yield a probability for each tt event to have 0, 
1, or > 2 6-tagged jets Q. 

We split the selected £+jets sample into subsamples 
according to the lepton flavor (e or //), jet multiplicity 
(3 or > 4 jets) and number of identified 6-jets (0, 1 or 

> 2), thus obtaining 12 disjoint data sets. We fit si- 
multaneously R and <j t i to the observed number of 1 b 
tag and > 2 b tag events, and, in b tag events with 

> 4 jets, to the shape of a discriminant T> that exploits 
kinematic differences between the background and the tt 
signal and which is described in detail below. As the 
signal-to-background ratio is about five times smaller in 
events with b tags and three jets we do not use a dis- 
criminant for that subsample. 

The dominant background is the production of W 
bosons in association with heavy and light flavor jets 
(W+jets). Smaller contributions arise from Z+jets, di- 
boson and single top quark production. Multijet events 
enter the selected sample if a jet is misidentified as 
an electron (e+jets), or a muon coming from either a 
semileptonic heavy quark decay or an in-flight pion or 
kaon decay in a light flavor jet appears isolated (/z+jets). 

We model the signal and backgrounds other than mul- 
tijet using a Monte Carlo (MC) simulation. The pro- 
cesses W+.jets and Z+jets are generated with the ALP- 
GEN 2.05 [12|] leading-order event generator for the multi- 
parton matrix element calculation and PYTHIA 6.323 [HI 
for subsequent parton showering and hadronization. Di- 
boson samples are generated with pythia and single 
top quark production is modeled using the SINGLETOP 
[lj] generator. The tt signal is simulated with pythia 
for a top quark mass of mt op = 175 GeV and includes 
three decay modes tt — > W + bW~b, tt — > W + bW~q~i (or 
tt — > W + qiW~b) and tt — ► W + qiW~q~i, where qi denotes 
a light down-type (d or s) quark. These three decay 
modes are referred to as bb, bqi and qiqi. The gener- 
ated events are processed through the GEANT-based [HI 
simulation of the DO detector. The same reconstruction 



5 



TABLE I: Sample composition for the measured a t i and R - 
1. The uncertainties are statistical only. 



7Vj ots sample 



b tags 



1 b tag > 2 b tags 



W+]ets 
Multijet 
Other 

ti 

Total 
Observed 



1394.4 ± 48.4 
287.4 ± 35.9 
254.0 ± 1.9 
109.7 ± 0.4 

2045.5 ± 60.3 
2050 



102.5 ± 
28.1 ± 
29.4 ± 
143.3 ± 
303.3 ± 
294 



3.5 
3.5 
0.2 
0.5 
5.0 



8.3 ± 
3.3 ± 
5.2 ± 
54.3 ± 
71.2 ± 
76 



0.3 
0.4 
0.4 
0.2 
0.5 



W+jets 
Multijet 
Other 
ti 

Total 
Observed 



188.2 ± 
66.9 ± 
62.2 ± 
83.8 ± 

401.1 ± 
389 



15.0 

9.9 

0.9 

0.3 

18.0 



17.3 ± 
6.6 ± 
8.0 ± 
126.4 ± 
158.3 ± 
179 



1.4 
1.0 
0.1 

0.5 
1.8 



1.8 ± 
0.8 ± 
1.7 ± 
64.2 ± 
69.5 ± 
58 



0.1 
0.1 
0.0 
0.2 
0.3 



algorithm as for data is used. Additional corrections [10( 
are applied to the reconstructed objects to improve the 
agreement between data and simulation. 

The determination of the background composition 
starts with the evaluation of the multijet background 
for each jet multiplicity and lepton flavor before &-jet 
tagging by counting events in the corresponding control 
data samples and applying the matrix method We 
estimate the number of events with a lepton originating 
from a W or Z boson decay by subtracting the estimated 
multijet background from the observed event yield be- 
fore ^-tagging. We further subtract diboson, single top 
quark and Z+jets contributions, normalized to the next- 
to-leading order cross sections [16j. The remaining data 
events are assumed to come from tt and W+jets back- 
ground. In the first step of the fitting procedure used 
to extract a t i and R,we assume the tt contribution pre- 
dicted by the SM |17j. In every subsequent step, we it- 
eratively re-determine the expected number of tt events 
and re-evaluate the W+jets background. 

Since the probability to tag a tt event depends on 
the flavor of the jets, it depends on R. We estimate the 
acceptance and tagging probabilities for each of the three 
tt decay modes bb, bqi and qiqi. Thus, the probability for 
a tt event to pass our selection criteria and to have n 
6-tagged jets is: 



Ptotal 



R 2 A(bb)P t n (bb) + 2R(1 - R)A(bq l )Pl l (bq l ) 
(l-R) 2 A( qm )P t n ( qm ), 



where A (P™) describes the acceptance (tagging proba- 
bility) for each decay mode. Figure[lja) shows P t n as a 
function of R for tt events with > 4 jets and 0, 1 and > 2 
b tags. Table U presents the composition of the selected 
sample for the measured a t i and R = 1. 

The topological discriminant D [l(| exploits the kine- 
matic differences between tt and W+jets events to 
achieve a better constraint on the number of tt events 
in the subsample with > 4 jets and b tags. We select 



variables well-described by the background model in sam- 
ples of events with one or two jets that provide a good 
separation between signal and VF+jets background. The 
optimal set of variables is chosen to minimize the ex- 
pected statistical uncertainty on the fitted fraction of tt 
events. Due to the differences in acceptance and sample 
composition, the discriminants are constructed from dif- 
ferent sets of variables in the e+jets and /i+jets channels. 
In the e+jets channel we use five variables: the leading 
jet px, the maximum A1Z fiol ] between two of the four 
leading jets, A, Cm, an d I'm [HI- In the ^+jets chan- 
nel the discriminant is built from six variables: A, T>m, 
the scalar sum of the pr of the four leading jets and the 
muon, the scalar sum of the pr of the third and fourth 
jet in the event, the transverse mass of the vector sum 
of all jets, and the ratio of the mass of the three leading 
jets to the mass of the event, defined as the invariant 
mass of the vector sum of all four jets, the lepton from 
the W decay and the missing transverse energy coming 
from the neutrino. The sensitivity to soft radiation and 
to the underlying event is reduced by using only the four 
highest-pr jets for the kinematic variables. 

The discriminant function is built using simulated 
VF+jets and ti events. We evaluate it for each physics 
process considered in the analysis and build correspond- 
ing template distributions consisting of ten bins. For ti 
we obtain a distribution for each of the three decay modes 
(bb, bqi and qiqi). The shapes of the discriminant distri- 
butions for Z+jets, diboson and single top backgrounds 
are found to be similar to that of the VF+jets events and 
we use the latter to model them. We use a sample of 
data events selected by requiring that the lepton fail the 
isolation criteria to obtain the discriminant shape for the 
multijet background. 

We define a likelihood function as the product of Pois- 
son probabilities over all 30 subsamples and bins of the 
discriminant, where in each subsample the expected num- 
ber of events is estimated as a function of R and a t i- 
We include 12 additional Poisson terms to constrain the 
multijet background prediction in each subsample. The 
systematic uncertainties are incorporated in the fit using 
nuisance parameters [l||, each represented by a Gaus- 
sian term in the binned likelihood. In this approach, 
each source of systematic uncertainty is allowed to affect 
the central value of R and a t i during the likelihood max- 
imization procedure, yielding a combined statistical and 
systematic uncertainty. 

The result of the maximum likelihood fit is: 



R 

on 



0.97±o;^ (stat+syst) and 

(stat+syst) ± 0.50 (lumi) pb 



— O 1 n+0.90 

- O-LO-0.84 



for a top quark mass of 175 GeV. Figure[TJb,c) compares 
the distribution of the data to the sum of predicted back- 
ground and measured tt signal for R = 0.97. We observe 
no significant dependence of R on mt op within ±10 GeV 




FIG. 1: (a) Fractions of events with 0, 1 and > 2 b tags as a function of R for ti events with > 4 jets; (b) predicted and 
observed number of events in the 0, 1 and > 2 b tag samples for the measured R and a t i for events with > 4 jets and (c) 
predicted and observed discriminant distribution in the b tag sample with > 4 jets. 



around the assumed value while a t i changes by +0.09 pb 
per 1 GeV within the same range. We find a correlation 
between R and a t i of -58%. Table HT1 summarizes the sta- 
tistical and leading systematic uncertainties on R and a t i 
excluding the 6.1% uncertainty on the integrated lumi- 
nosity [20( ■ The contribution of each individual source of 
uncertainty is estimated by fixing all but the correspond- 
ing Gaussian term in the fit. The statistical uncertainty 
is obtained from the fit with all Gaussian terms fixed. 

The total uncertainty on R is about 9%, compared 
to 17% achieved in the previous measurement Q. The 
largest uncertainty comes from the limited statistics. 
Since the 6-tagging efficiency drives the distribution of 
the events among the 0, 1 and 2 6-tag subsamples and 
is strongly anti-correlated with R, the systematic uncer- 
tainty is dominated by the 6-tagging efficiency estima- 
tion, responsible for ^90% of the total systematic uncer- 
tainty. 

The total uncertainty on a t t, excluding luminosity, is 
~10.5%, representing a 30% improvement over the pre- 
vious measurement [7] performed under the assumption 
of R = 1. In the latter, the primary 4.7% relative un- 
certainty comes from the 6-tagging efficiency estimation 
while in the current measurement it is reduced to 1.2% 
because a t i is much less sensitive to the variations of the 
6-tagging efficiency than R, and the two-dimensional fit 
takes advantage of this feature. 

We extract a limit on R and \Vtb\ following the 
Feldman-Cousins procedure [21[. We generate pseudo- 
experiments with all systematic uncertainties included 
for various input values of R (i?tmc) and apply the 
likelihood-ratio ordering principle. We obtain R > 0.88 
at 68% C.L. and R > 0.79 at 95% C.L., illustrated in 
Fig. [2] From R we determine the ratio of \Vtb\ 2 to the 



TABLE II: Summary of uncertainties on a t i and R. 



off-diagonal matrix elements to be 



\v tb 



> 3.8 at 



\Vts\ 2 + \Vtd\ 2 

95% C.L. Assuming a unitary CKM matrix with three 
fermion generations we derive |Vt&| > 0.89 at 95% C.L. 



Source 




(Pb) 


AR 




Statistical 


+0.67 


-0.64 +0.067 - 


0.065 


Lepton identification 


+0.32 


-0.27 


n/a 




Jet energy scale 


+0.32 


-0.23 


n/a 




TF+jets background 


+0.21 


-0.23 


n/a 




Mufti jet background 


+0.17 


-0.17 +0.016 - 


0.016 


Signal modeling 


+0.12 


-0.25 


n/a 




6-tagging efficiency 


+0.10 


-0.09 +0.059 - 


0.047 


Other 


+0.24 


-0.13 +0.015 - 


0.014 


Total uncertainty 


+0.90 


-0.84 +0.092 - 


0.083 



In summary, we have performed a simultaneous mea- 
surement of the ratio of branching fractions R and o t i 
yielding the most precise measurements to date, R = 
0.97±2;81 (stat+syst) and a t{ = 8.18±g;|2 (stat+syst) ± 



0.50 (lumi) pb, both in good agreement with the SM. 
This measurement of R will be a key ingredient in a fu- 



ture model-independent direct determination of the \V t 



i<i\ 




FIG. 2: The 68% (inner band), 95% (middle band) and 99% 
(outer band) C.L. bands for -Rtruc as a function of R. The 
dotted black line indicates the measured value R = 0.97. 
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CKM matrix elements. 
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